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We study superflow decay via quantum phase slips in trapped one-dimcnsional (ID) quantum 
gases through dipole oscillations induced by sudden displacement of the trapping potential. We find 
the relation between the damping rate of the dipole oscillation G and the phase-slip nucleation rate 
r as G oc r/u, where v is the flow velocity. This relation allows us to show that damping of ID 
Bose gases in optical lattices, which has been extensively studied in experiment, is due to quantum 
phase slips. It is also found that the damping rate versus the flow velocity obeys the scaling formula 
for an impurity potential even in the absence of an explicit impurity. We suggest that the damping 
rate at a finite temperature exhibits a universal crossover behavior upon changing the flow velocity. 
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Systems of optical lattices loaded with ultracold gases 
have offered unique opportunities for the studies of cor- 
related many-body physics in one dimension (ID), owing 
to their extraordinary controllability and cleanness [l|, ■ 
In typical experiments, one creates an array of ID gases 
by focusing a strong 2D optical lattice to a 3D gas, and 
the transverse confinement is widely controllable so that 
thermal and quantum motions in the transverse direction 
are completely frozen. With such extremely ID quan- 
tum gases, recent experiments have revealed intriguing 
phenomena that are in stark contrast to higher dimen- 
sions, including Tonks- Girardeau gases 0] and their 
non-ergodic dynamics Q , a possible fermionic superfluid 
of the Fuldc-Fcrrell-Larkin-Ovchinnikov type Q, the pin- 
ning Mott transition by a shallow optical lattice , and 
strong suppression of superfluid transport 04lll- 

Transport of trapped ultracol d g ases through peri- 
"T^l, single-barrier or random potcn- 



odic 

tials [111 U^l has been often investigated by suddenly dis- 
placing a parabolic trap to induce a dipole oscillation 
(DO) and observing its damping. As for the transport of 
ID Bose gases, it has been found that the DO in the pres- 



ence of an optical lattice or random potential 11 
the axial direction is significantly damped even in the su- 



171 have 



perfluid state. Previous theoretical studies 
suggested that this apparent contradiction, namely dis- 
sipative flow in a superfluid, can be interpreted as a con- 
sequence of phase slips (PS), in which thermal or quan- 
tum fluctuations cause the phase of the superfluid order 
parameter to unwind, leading to the dissipation of flow. 
This interpretation, if affirmative, could provide a unified 
view for superfiuidity in ID, given that the concept of PS 
is central also to the understanding of ID superfiuidity 
and superconductivity in several other condensed-matter 
systems, such as liquid '^He in ID nanopores [13, flit - 
metallic nanowires |2C 



2l| , and single- walled carbon nan- 



otubes [22-24|. Moreover, thanks to the fiexible control- 
lability of optical lattice systems, it would open up new 



possibilities for more thorough studies of PS. However, 
relating explicitly the damping of DO to PS is highly 
non-trivial and has never been done thus far, because of 
difficulty in analyzing PS under a non-uniform trapping 
potential. Indeed, despite a number of previous numeri- 
cal works on damped DO of ID gases [251434 1. the inter- 
pretation in terms of PS has never been corroborated. 

In this Letter, we study the DO dynamics of trapped 
ID superfluids in connection with PS. Through quali- 
tative consideration on energy loss during the damping 
and exact numerical simulations with time-evolving block 
decimation (TEBD) method [HI at zero temperature, we 
find a parameter region where the damping rate G and 
the PS nucleation rate F satisfy the following simple re- 
lation as a function of the flow velocity u, 



G{v) (X r{v)/v. 



(1) 



We emphasize that since the damping rate is a major ex- 
perimental observable [3,[^,[ll|, the relation ([1]) allows for 
analyzing the PS nucleation rate in experiment. Using 
this relation and TEBD, we show that the damping rate 
in ID Bose gases in an optical lattice obeys the power-law 
formula derived from the nucleation rate of a quantum 
PS (QPS), thus numerically confirming the PS scenario. 
The exponent of the power-law is found to coincide with 
that for an impuri ty p otential (36l- 38| rather than for a 



periodic potential [l7[, although there is no explicit im- 
purity. We also discuss the effects of finite temperatures 
to suggest a universal behavior of the damping rate in ID 
superfluids that can be considered as a single-component 
Tomonaga-Luttinger (TL) liquid. 

We describe ID Bose gases of the total particle number 
N by means of the following ID Bose-Hubbard model, 



H = -J^(6]6,+i+h.c.) + -^n,(n,-l) 

j 3 

+y ^ + ^ - Xjdf + A<5,- o] h,,{2) 

j 



2 




FIG. 1: Numerical data for the Bose-Hubbard model Eq. ([2]) 
in the hardcore limit {U — i> oo). We set A*' = 31, V/J = 
-1.4, n/J = 0.00032, A/J = 1, and xo/d = 8. (a): The 
solid, dashed, and dashed-dotted lines represent the density 
distribution rij = {hj), the external potential at t < 0, and 
that at t > 0. (b) and (c): The time evolution of the center 
of mass position Xc.m.(t) and velocity Wc.m.(i). 



where bj denotes the annihilation operator on the jth site 
and fij = bjbj. J, U, and V are the hopping energy, the 
onsite interaction, and the nearest-neighbor interaction, 
respectively. The last term in Eq. ^ means the external 
potential that consists of a parabolic trapping potential 
and a single impurity, where fi is the trap curvature, Xc 
the position of the trap center, d the lattice spacing, and 
A the impurity strength. While the system can be in an 
insulating state for certain values of the parameters, we 
hereafter assume the superfluid state because our main 
interest is in superflow decay via PS. 

We use the TEBD method (s^l to calculate the ground 
state and the exact quantum dynamics of the DO of 
Eq. ([2]). TEBD is a quasi-exact numerical method for 
computing time evolution of a many-body wave function 
of ID quantum lattice systems, and it can prec isely de- 
scribe quantum fluctuations causing QPS |l7l . Isol . lioj . 
To prepare an initial state for the real-time evolution, 
we first compute the ground state of Eq. ^ with the 
trap displacement Xc = xq (> 0) via imaginary time 
evolution. In Fig. [ija), we show a typical external po- 
tential (dashed line) and density profile (solid line) of 
the initial state. With this initial state, we displace 
the trap to = at t = as shown by the dashed- 
dotted line in Fig. [Ifa) and compute the real-time evo- 
lution. In Figs. [Ijb) and (c), we show an example of 
the time evolution of the center of mass (cm.) position 
Xc.m. = N~''-dJ2j ji^j) and velocity Vc.m. = ic.m., which 
exhibit a damped DO. We extract the damping rates us- 
ing the formula, G = In {Aq /Ai ) /ti . As indicated in 
Fig. [IJb), Aq and Ai arc the initial amplitude of Xc.m./d 
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FIG. 2: Damping rates G versus the maximum flow velocity 
I'max for the hardcore Bose-Hubbard model with A'' = 31 and 
several values of A/J. We set {V/J,n/J) = (0,0.001) (a) and 
(-1,0.0005) (b). The thin solid lines represent l/(4ii) (in 
the unit of J/fi) at A = and xo = d. The thick solid lines 
represent lOGo. The dashed lines represent Umax at xo = d. 
The dotted lines represent Vc/5. Numerical fltting with a 
function G{v) — Cv^ is made for the data within the shaded 
region in order to extract the prefactor C and the exponent rj, 
and each dashed-dotted line is best fit to data for each A/J. 

and the amplitude after the half period ti. 

Before numerically verifying the relation ([1} between 
G and F, let us present a qualitative explanation that 
provides intuitive understanding of the relation. For this 
purpose, we express the oscillation-energy loss through 
the damping of the first half period in the following two 
ways. The first one is in terms of the lost potential en- 
ergy, -Eioss = ^-^'^^(^0 ~ ^i): where oj is the oscillation 
frequency and M the total mass. Assuming that the 
damping is so weak that i5 = 1 — Ai/Aq ^ 1, or equiva- 
Icntly Gti <C 1, the energy loss can be rewritten as 

^loss =^ MiujAofS ~ Mvl^^Gh . (3) 

Note that we used S ~ Gti and Wmax — i-^Aq to derive 
Eq. dSl), where "^max is the maximum cm. velocity as indi- 
cated in Fig.[ljc). Thus, E\oss is expressed with G. The 
other way is to use the Joule heat, -Eioss = P ^ ti, where 
P = RI^ is the power, R the resistance, / ~ J^avo^^max 
the particle current, and riavc the average density. As- 
suming that the main source of the resistance is due to 
PS, the resistance can be related to the nucleation rate 
as i? = 2-KhT / 1 [4l|. Using this relation, one obtains 

E\oss ^ STrftnavcWmaxPtl, (4) 

which connects the energy loss and the nucleation rate. 
Equating the right-hand side of Eq. (jSj with that of 
Eq. ([U leads to the relation, G ^ 2nhn^yc/^I x P/wmax, 
which agrees with the relation ([1]). 

Since the above explanation of the relation ^ is 
only qualitative, more quantitative corroboration is de- 
manded. Hence, we analyze the DO in the hardcore bo- 
son limit {U — > oo) using TEBD, in order to provide ac- 
curate numerical verification of the relation ([TJ . In this 
limit, as long as —2 < V/J < 0, low energy physics of 



3 



the system is known to be well described by the TL liq- 
uid model [i^. Previous analytical studies have shown 
that the nucleation rate of a QPS in the TL liquid with a 
single impurity exhibits the following power-law behav- 



ior with respect to v as Fgi 



,2if-l 



V <ti Vc [3^ 37 1 . Here Vc is the mean- field critical veloc- 
ity and K is the TL parameter the inverse of which, 
namely 1/K, quantifies the strength of quantum fluctu- 
ations. To hold K under control, we fix iV = 3f , and ad- 
just fl/ J depending on V/ J such that rimax — 0.5, where 
nmax is the maximum density [see Fig. (Ha)]. In such 
a situation, the analytical expression at half filling Q, 
K = 7r/[27r — 2 arccos (57)], is approximately valid, and 
K can be controlled by changing only V/J. 

If the relation ([T]) is correct, the damping rate should 
obey G cx v^^^"^. To corroborate this, we plot in 
Figs. [2ja) and (b) G versus Umax for V/ J = Q and —I, 
taking different values of A. Note that we vary xq to 
control Wmax. As indicated by the shaded area in Fig. [21 
we find the parameter region in which the damping rate 
safely obeys the power-law formula. This region is de- 
termined by the following four conditions: i) Gti < 1/4, 
ii) G > lOGo, where Gq is the damping rate at A = 0, 
in) a;o > d, and iv) Umax < Vc/5 [i^l- We recall that the 
relation ^ is supposed to be valid when Gti <C 1 and 
the source of the damping is mainly due to PS. While 
the first condition obviously corresponds to the former 
requirement, ii) and iii) stem from the latter. As for ii), 
we see from the black circles in Fig. [2] that there is small 
but finite damping because of dcphasing effects even in 
the absence of an impurity that induces QPS [l^]. To 
distinguish the QPS from the dephasing, the damping 
at A > has to be much larger than that at A = 0, 
thus requiring ii). The condition iii) is necessary be- 
cause otherwise a mismatch of the initial density and 
the displaced trap causes additional damping or revival 
of the DI that blurs the QPS effects. The last condition 
has to be satisfied to validate F c>c v'^^~^, as mentioned 
above. These four conditions are indicated by the thin- 
solid, thick-solid, dashed, and dotted lines in Fig. [2l 

The damping rates in the shaded area surrounded by 
the four lines indeed exhibit the power-law behavior. By 
fitting a function G{v) = Cv^ to the data in the area, 
we extract the exponent 77 and the prefactor G, where 
G = KG/ J and v = hvmax/iJd). In Figs. |3|a) and 
(b), we plot T] and C ioi V/J = and —1 as functions 
of A/ J. When A/ J increases, rj is almost constant and 
nearly equal to 2K — 2, and G quadratically increases for 
A < J. This is consistent with the previous results that 



for any A when 



..2K-1 



holds for any A 



and that Fs 



A^ 



for small A 136|. In Figs. [31Jc), we plot 77 versus V/J for 
different values of A/ J. There we see that the exponents 
agree very well with the expected value, i.e., 77 = 2 A' — 2 
that is represented by the solid line in Fig. |3ljc) . 

Having corroborated the relation ([T]) in both qualita- 
tive and quantitative manners, we now consider the case 
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FIG. 3: The exponent ?7 and the prefactor C for the hardcore 
limit with A'^ = 31. (a) and (b): 77 and C as functions of 
X/J for {V/J, ft/ J) = (0,0.001) (red circles) and (-1,0.0005) 
(black squares). The red dashed and black solid lines in (a) 
represent rj = 2K — 2 for V/ J — and —1. The red dashed 
and black solid lines in (b) represent a parabolic function 
/(A/J) = a{X/jf for V/J = and -1 with the constant 
a determined such that the lines pass on the data points at 
A/J = 0.6. (c): 77 versus V/J for several values of A/J. The 
solid line represents rj = 2K — 2. 
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FIG. 4: The case of the softcore Bose-Hubbard model {U < 
00) with V = and A = 0. (a): G versus TJmax for U/J — 3.2 
and fi/J = 1/900. The thin solid, thick soUd, dashed, and 
dotted hnes represent l/(4ti) at xo = d, lOGiow, ^max at 
xo = d, and Vmax = Vc/8. The dashed dotted hne represents 
the best fit to the data inside the shaded region, (b): 77 versus 
U/J. The solid and dashed lines represent 2K — 2 and 2K — 
3. The shaded regions mean the error bars of these lines 
originating from the errors in numerically evaluating K. 



of softcore bosons {U < 00) without the nearest neigh- 
bor interaction and the impurity. This case is of direct 
relevance to the experiments of Refs. 



11 



where 

the damped DO of ID Bose gases in optical lattices has 
been studied, and is of great importance for understand- 
ing whether the damping observed in the experiments is 
due to QPS. To address QPS effects, we choose TV and ft 
such that 1 < ?7.max < 2. In this situation, there exist the 
regions of nj ~ 1, where the underlying lattice structure 
induces strong umklapp scattering leading to QPS. Since 
the nucleation rate of such a QPS obeys Fprd 



OC V 



2K-2 



when V < Vc [17|, we naively speculate G cx v"^^'^. We 
take U < Uc for the system to be safely in the snperfluid 
state, where Uc — 3.3J [3l is the Mott transition point. 
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FIG. 5: (a): Schematic crossover diagram for the damping 
rate G as a function of the temperature T. (b): Sketch of the 
universal behavior in G versus w at a finite temperature. 



In Fig. mja), wc plot G versus 



for N 



99, 



n/J = 1/900, and U/ J — 3.2. Wc again find the param- 
eter region in which the damping rate exhibits a power- 
law behavior. While the conditions i) and iii) remain the 
same as the hardcore boson case, ii) and iv) are modified 
as ii') G > lOGiow and iv') Wmax < Vc/8, where Glow is 
the damping rate for rimax < Ij C-g-i = 37 is taken to 
be compared with the case shown in Fig.^Ja). The condi- 
tion ii') guarantees that the main source of the damping 
is not the dephasing but the PS occurring at rij ~ 1. The 
reason for the latter modification will be discussed later. 

By fitting a function G{v) = Gv^ to the data in the 
shaded area, we extract r] for {N,n/J) = (99,1/900) 
and (67,1/400), and plot them as functions of U/J in 
Fig. m^b). rj agrees with the value expected for an impu- 
rity potential, 2K — 2 (solid line), rather than that for a 
periodic potential, 2K — 3 (dashed line) This dis- 

agrees with the naive speculation mentioned above and 
is counter-intuitive in the sense that there is no explicit 
impurity in the system. However, it can be interpreted 
as follows. When a Bose gas is trapped in a combined 
parabolic and periodic potential and rimax > 1, there 
are two narrow regions where rij ~ 1. In these regions, 
the umklapp process is the most relevant, and the trans- 
port is strongly suppressed. Hence, the unit-filling re- 
gions move slower than the other parts of the gas, and 
act as impurities for the other parts. This interpretation 
also explains the requirement iv') in the sense that the 
presence of the effective impurities may lower Vc- 

On the basis of the finding that the damping rate is 
related to the PS nucleation rate for an impurity po- 
tential both in the presence and the absence of an ex- 
plicit impurity, we suggest that there are the follow- 
ing three distinct regimes with respect to the damping 
due to PS at finite temperatures, which are illustrated 
in Fig. EKa). (A) When the temperature is as low as 



ksT <C Ejv/vc, the PS is caused by pure quantum tun- 
neling and the damping rate obeys G oc v^^~^, as dis- 
cussed above. Here Ej = hu/{\^d) is the Josephson 
plasma energy and u is the sound velocity. (B) In the 
intermediate temperatures, Ejv/vc <S; <g; SF, the 
PS occurs due to the thermally assisted quantum tun- 
neling j47| and F oc vT -2 [ii, [13, corresponding to 
G cx T'^K~^ , where 5F is the free energy barrier separat- 
ing two neighboring winding-number states. (C) When 
UbT ^ 5F, the thermal activation process becomes dom- 
inant [m and G cx e-^^/(^^'^\ 

As sketched in Fig. [^b), when k^T <t: SF ^ Ej, the 
crossover between the regimes (A) and (B) can be in- 
duced by changing the flow velocity with a fixed temper- 
ature. The two regimes are separated by the crossover 
velocity that is ^ Vc x kBT/Ej. Given that one can 
achieve Ej/kB - 30nK [i| and T - 4nK ^ in cur- 
rent experiments, it is likely that the crossover can be 
observed. Since the main feature of this crossover is de- 
termined only by the TL parameter K and the sound 
velocity u, we conjecture that it can be applied univer- 
sally to the damping of DO in ID quantum gases that 
can be effectively described as a single-component TL liq- 
uid. Examples include not only the two cases addressed 
above, but also paired or counterflow superfluid states of 



34| gases. 



two-component Bose [33[ or Fermi [3^ 

In conclusion, we have studied damped DO of trapped 
ID quantum gases in the presence of an impurity po- 
tential or an optical lattice from a perspective of PS. We 
have connected the damping of DO and the PS nucleation 
through the relation ([T]). Combining this relation with 
the TEBD simulations of the ID softcore Bose-Hubbard 
model, we found that in certain parameter regions the 
damping rate algebraically grows with increasing the flow 
velocity as expected from the nucleation rate of a QPS. 
This result strongly supports the interpretation that the 
strong suppression of superfluid transport observed in the 
experiments . [in is due to QPS. We also suggested a 
universal damping behavior at finite temperatures, which 
can be tested in future experiments. 
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manuscript and valuable comments. Useful discussions 
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